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We compute in detail the rate at which energy injected by dark matter annihilation heats and ion- 
izes the photon-baryon plasma at 2 ~ 1000, and provide accurate fitting functions over the relevant 
redshift range for a broad array of annihilation channels and DM masses. The resulting pertur- 
bations to the ionization history can be constrained by measurements of the CMB temperature 
and polarization angular power spectra. We show that models which fit recently measured excesses 
in 10-1000 GeV electron and positron cosmic rays are already close to the 95% confidence limits 
from WMAP. The recently launched Planck satellite will be capable of ruling out a wide range of 
DM explanations for these excesses. In models of dark matter with Sommerfeld-enhanced annihi- 
lation, where {av) rises with decreasing WIMP velocity until some saturation point, the WMAP5 
constraints imply that the enhancement must be close to saturation in the neighborhood of the 
Earth. 

PACS numbers: 95.35.-|-d,98.80.Es 
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I. INTRODUCTION 

Dark matter (DM) annihilation around the redshift 
of last scattering {z ^ 1000) can modify the observed 
temperature and polarization fluctuations of the CMB, 
which have been measured to high precision by exper- 
iments such as WMAP, ACBAR and BOOMERANG 
[I) S Hi ■ Other indirect astrophysical probes of DM an- 
nihilation must contend with the complexities of Galac- 
tic astrophysics - for example, the DM distribution and 
dumpiness, ISM density, magnetic field strength and de- 
gree of tangling, Galactic photon energy density, etc. All 
of these are complex processes with significant uncertain- 
ties; in contrast, the mechanisms by which DM annihi- 
lation modifies the CMB are relatively simple and well 
understood. 

DM annihilation injects high-energy particles into the 
IGM [7l| , which heat and ionize neutral hydrogen as they 
cool. This ionizing energy does not generally change 
the redshift of recombination, but does alter the resid- 
ual ionization after recombination. The increased ion- 
ization fraction leads to a broadening of the last scat- 
tering surface, attenuating correlations between temper- 
ature fiuctuations. The low-£ correlations between polar- 
ization fiuctuations, on the other hand, are enhanced by 
the thicker scattering surface. 

These effects of WIMP annihilations on recombina- 
tion have been studied previously, with significant effects 
on WMAP for (cta") ~ (IQ-^e//) x (AfoM/llGeV/c^)) 
cm^/s and on Planck for (ctav) ^ (10"^''//) x 
(MDM/(lGeV/c^)) cm'^/s, where / is an unknown pa- 
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rameter of order - 0.01 -14:]. For O(TeV) WIMPs, 
these cross sections were considerably larger than the 
standard thermal cross sections expected at the time, but 
in light of recent cosmic-ray experiments which motivate 
larger annihilation cross sections, it is important to re- 
consider these limits and quantify the free parameter /. 

The PAMELA, ATIC, PPB-BETS, Fermi and H.E.S.S 
experiments have observed unexpected features in the 
electron and positron cosmic-ray spectra, at energies of 
10-1000 GeV. PAMELA has measured the positron 
flux ratio (j){e'^) / {(j){e'^) + (/'(e~)) and found a sharp rise 
starting at 10 GeV and continuing up to ~ 100 GeV. This 
confirms previous indications of an excess by HEAT Q 
and AMS-01 0, but at much higher confidence. The 
ATIC balloon experiment Q has measured the spectrum 
of e+ + (ATIC cannot distinguish positrons from elec- 
trons) from 20-2000 GeV, and finds a broad excess at 
300 — 800 GeV, in agreement with the similar excess ob- 
served by PPB-BETS ^. The Fermi ^ and H.E.S.S 
[m experiments have measured a similar but somewhat 
smaller e"*" -I- e~ excess in the 300 — 1000 GeV energy 
range, and have not confirmed the peak and sharp cutoff 
observed by ATIC around 700 GeV. However, in combi- 
nation, these results argue for a new primary source of 
high-energy electrons and positrons. 

The annihilation of WIMPs in the Galactic halo has 
been proposed as this new source, with the other major 
contender being e~^e~ pairs from pulsars fl^. ITsl. [l^ [lq| . 
While weak-scale dark matter annihilating in the Galac- 
tic halo can generically produce annihilation products 
with energies in the hundreds of GeV, including elec- 
trons and positrons, most DM annihilation channels give 
rise to insufficiently hard electron spectra. For ener- 
gies up to 100 GcV, PAMELA observed no antiproton 
excess [l6j . strongly constraining annihilation channels 
that produce copious antiprotons. Gamma ray bounds 
also constrain the production of neutral pions (either di- 
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rectly or as products of a hadronic cascade). Taking these 
constraints into account, for dark matter masses in the 
10 GeV - 1 TeV range, the WIMP must annihilate pri- 
marily to leptons in order to fit PAMELA observations 

Even for these primarily leptonic annihilation chan- 
nels, the annihilation cross section required to fit the 
observed excesses is significantly higher than the ther- 
mal rehc cross section 0, |20|. For a WIMP with 
^ TeV mass, as suggested by the measured charge- 
undifferentiated electron spectrum, a boost factor of 2-3 
orders of magnitude is required, depending on the an- 
nihilation channel. Lower masses allow smaller boost 
factors, but a 100 GcV WIMP still requires a boost 
factor of ~ 10 or higher relative to the thermal relic 
cross section, to produce the PAMELA positron excess 
(and of course, cannot produce the observed ATIC/PPB- 
BETS/Fermi/H.E.S.S excesses). 

Several authors have suggested models of thermal relic 
DM where {<jav) rises with decreasing WIMP velocity 
as a consequence of a Sommerfeld enhancement, and 
thus the present-day DM annihilation cross section is 
considerably higher than the thermal relic cross section 
(aAv) ~ 3x 10-26 cmVs [H [H, IH, [H, [H . Sommerfeld 
enhancement occurs when some particle with mass much 
smaller than the WIMP mass mediates a long-range at- 
tractive force between the WIMPs [H, For heavy 
WIMPs (~ 10 TeV), these force carriers can simply be 
the Standard Model W and Z bosons [2l[. Other mod- 
els invoke new light dark-sector particles [H, [H, [2^ to 
mediate a Sommerfeld enhancement for weak-scale DM 
(~ 100 — 1000 GeV). In general, Sommerfeld-enhanced 
models suggest that a much higher annihilation cross sec- 
tion around the redshift of recombination may be possi- 
ble. 

Sommerfeld-enhanced dark matter models with weak- 
scale DM masses and GeV-scale light force carriers are 
especially attractive for explaining the observed cosmic 
ray anomalies [23l . [23 |. If the DM annihilates to the 
light force carriers, then decays into protons/antiprotons 
(and gauge bosons, etc) are kinematically forbidden, and 
the annihilation products are predominantly hard lep- 
tons [l!]. The electron and positron spectra produced 
by such annihilation channels provide excellent fits to 
the ATIC and PAMELA e+e" spectra, up to a boost 
factor attributed to the Sommerfeld enhancement [20| . 
and similarly fit the Fermi data well (23j- Such mod- 
els have also been proposed to explain anomalous results 
from the DAMA [s^ experiment with inelastic WIMP- 
nucleon scattering [31,^, 33.], and the INTEGRAL [34] 
signal with inelastic WIMP- WIMP scattering [H, [3^ . 
which naturally arises in this framework [23|. We adopt 
the terminology of [35] , and refer to these models as "ex- 
citing dark matter" (XDM). CMB probes of DM annihi- 
lation are particularly well suited to constraining mod- 
els with Sommerfeld-enhanced annihilation, since the 1 /v 
scaling of (aAv) results in a large annihilation cross sec- 
tion when the DM is diffuse and cold. 



In order to constrain specific models of DM annihi- 
lation, it is essential to understand the efficiency with 
which energy from DM annihilation heats and ionizes 
the IGM, as it is this deposited energy which perturbs 
the ionization history and hence the CMB. Several pre- 
vious studies have investigated the effect of dark matter 
annihilation/decay on the ionization history at recom- 
bination and during the cosmic "dark ages" [4, 37, 3J, 
m, EqI, SH, US]- However, these works have generally ei- 
ther assumed that the energy from annihilation/decay 
is deposited promptly with some redshift-independent 
efficiency, leaving this efficiency factor / as a free pa- 
rameter (degenerate with varying the annihilation cross 
section), or alternatively characterized the universe as 
opaque except for photons with energies lying in some 
"transparency window" . 

We calculate in detail the rate at which energy from 
dark matter annihilation is deposited into the IGM, via 
interaction of the annihilation products with the photon- 
baryon plasma. This allows us to apply previously de- 
rived constraints on f{aAv) to the DM model itself; in 
particular, we can directly apply the WMAP5 constraints 
to models proposed to explain the observed cosmic-ray 
anomalies. 

A similar calculation has been performed by [43|, and 
employed in [3] to discuss the effect of light dark mat- 
ter decays on the 21cm line. However, this analysis 
applied only to decay and annihilation of light (< 10 
MeV) dark matter, where high energy photon cooling 
mechanisms such as pair production and photon-photon 
scattering could be neglected, and also assumed a near- 
monoenergetic spectrum. In contrast, our calculation is 
valid for arbitrary spectra of the DM annihilation prod- 
ucts, for dark matter masses up to TeV scales. Our anal- 
ysis also improves on the prior calculation by taking into 
account the changing ionization fraction of the universe 
around recombination, and the energy injection from DM 
annihilation products at redshifts greater than z ~ 1100, 
which cool slowly and eventually deposit their energy at 
a lower redshift. 

Section |TT] details our numerical calculation of the en- 
ergy deposition from DM annihilation, while Section IIIII 
presents our results for the redshift-dcpendent efficiency 
factor f{z), for an array of annihilation channels. Sec- 
tion IIVI discusses WMAP5 constraints on the annihila- 
tion cross section for various annihilation channels, the 
implications for DM models with Sommerfeld-enhanced 
annihilation, and the prospect that Planck [45| can rule 
out dark matter annihilation as an explanation for the 
observed cosmic-ray anomalies. 



II. ENERGY LOSS PROCESSES FOR DM 
ANNIHILATION PRODUCTS 

Depending on the DM model, the dark matter may 
annihilate to a wide range of particles: gauge bosons, 
charged leptons, neutrinos, hadrons, or more exotic 
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states. These annihilation products may subsequently 
decay or interact with the IGM, producing showers of 
pairs, protons and antiprotons, photons and neutrinos. 
Neutrinos are stable and weakly interacting, so they es- 
cape and their energy is lost, while protons are highly 
penetrating and poor at transferring energy to the IGM 
[37| . Heating and ionization of the IGM occurs primarily 
through the electrons, positrons, and photons injected as 
a result of DM annihilation, either directly or by subse- 
quent decays and interactions of the primary annihila- 
tion products. Positrons behave identically to electrons 
at high energies, while at low energies they thermalize, 
form positronium and annihilate into photons. Thus the 
problem of computing the energy deposition efficiency 
reduces to calculating the evolution of the photon and 
electron spectra. 

The specific annihilation channels we focus on in this 
paper are motivated by the anomalous excesses observed 
by PAMELA and ATIC, and involve annihilation either 
directly to charged lepton pairs, or to light dark-sector 
states which decay to charged lepton pairs [l^. In the 
latter case, the same light dark-sector states mediate a 
Sommerfeld enhancement. We also compute the energy 
deposition efficiency for two benchmark masses (200 GeV 
and 1 TeV) for WIMPs annihilating via other SM chan- 
nels. 



A. Electron cooling and energy deposition 
mechanisms 

Energy loss processes for high energy electrons have 
been considered by a number of authors [13, [H, ■ The 
primary cooling mechanisms are inverse Compton scat- 
tering on CMB photons at high energy (7 ^ 1), and 
collisional heating, excitation and ionization at low en- 
ergy (see Appendix [Bl for the detailed cross sections). 

In all cases, electrons deposit their energy into the 
IGM, or produce high-energy photons via inverse Comp- 
ton scattering on the CMB, on timescales short compared 
to the Hubble time. For high-energy electrons, most of 
the energy goes into inverse Compton scattered photons, 
converting the injected electron spectrum to an effective 
injected photon spectrum. It is important to take into 
account repeated inverse Compton scatterings, as only 
in the extreme Klein-Nishina limit do electrons lose all 
their energy in the first scattering, so high-energy elec- 
trons can produce many low-energy photons via repeated 
scatterings. Once these electrons have cooled to low en- 
ergies where ICS is no longer the dominant process, they 
- like the electrons injected at low energy - rapidly de- 
posit their remaining energy to the IGM by ionization, 
excitation and heating. The positrons produced by DM 
annihilation behave identically to electrons at high ener- 
gies; at the lowest energies they form positronium and 
then annihilate, producing the usual positronium spec- 
trum with a contimmm up to 511 keV and a line at the 
cutoff. 



In order to determine the fraction of annihilation power 
deposited to the IGM as a function of redshift, we must 
compute the amount of energy deposited to the IGM and 
the spectrum of photons produced by these processes, 
when the energy in the injected electron spectrum has 
been completely depleted. The relative rates of the var- 
ious energy loss processes change with redshift, so we 
perform this calculation for an arbitrary injected elec- 
tron spectrum at each of the relevant redshifts. As all 
these processes occur on timescales much faster than the 
Hubble time, we do not consider redshifting effects, and 
so the fraction of energy loss into each mode can be de- 
termined simply by comparing their rates. For electrons 
of each energy, we record the fraction of energy that goes 
into created or upscattered photons and energy deposi- 
tion to the IGM, and the fraction retained by the electron 
as it downscatters to lower energies. 

We assume that below a certain kinetic energy 
(presently set to 250 eV) electrons efficiently deposit all 
their kinetic energy, cooling completely without produc- 
ing photons by ICS. We also make the approximation 
that after a single ionization, an electron deposits all its 
energy, since in the energy range where ionization domi- 
nates inverse Compton scattering, electrons tend to lose a 
large fraction of their energy to ionization, and the result- 
ing secon dary electrons efficiently deposit their energy to 
the IGM [U). 

For each initial electron energy, we can then determine 
the eventual partition of the electron's energy between 
produced photons and energy deposited to the IGM, 
when the electron has completely cooled and thermal- 
ized (and possibly annihilated). We use an inductive ap- 
proach, starting with the lowest energy bins and working 
up in energy. Since we have recorded the fate of all elec- 
trons with lower initial energy, after the primary electron 
has downscattered once we can employ the previously 
calculated cooling histories. 

Once the energy partitioning has been determined for 
all electron bins, integrating this result over the injected 
electron spectrum yields the prompt deposited energy 
and the photon spectrum produced by electrons injected 
at the given redshift. Determining the annihilation power 
that is eventually deposited to the IGM then becomes a 
question of how the photons deposit their energy. The 
dominant photon energy loss processes are not much 
faster than H{z) over the entire relevant energy range, 
so redshifting must be taken into account. 

Applying the same techniques to electrons and 
positrons from pair production, or electrons from Comp- 
ton scattering, allows us to view these mechanisms as 
converting a given photon spectrum to a new photon 
spectrum, plus some rapidly deposited energy. In this 
way we "integrate out" the fast electron cooling processes 
and reduce our numerical problem to computing the evo- 
lution of the photon spectrum only. 
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FIG. 1: A comparison of the photon coohng time to the Hub- 
ble time at z = 1000, for different photon energies. The 
dominant processes (in order of increasing energy) are ioniza- 
tion, Compton scattering, pair production on the H/He gas, 
photon-photon scattering, and pair production on the CMB. 
All the curves assume a He mass fraction of 1/4, with a density 
of 2.57 X 10~^ amu / cm^ today. The dotted curve shows pair 
production on a neutral IGM, the dashed curve shows pair 
production on a fully ionized IGM, and the dashed-dotted 
curve represents pair production on the CMB. This figure 
updates Fig. 1 in ,4], which had an error leading to cooling 
times approximately a factor of three longer. 



B. Energy deposition from photons 

The interaction of photons with the IGM was consid- 
ered in detail by 0, |4i|, who find that the dominant 
processes (ordered by increasing photon energy) are pho- 
toionization, Compton scattering, pair production off nu- 
clei and atoms, photon-photon scattering, and pair pro- 
duction off CMB photons. The cross sections and spectra 
for these processes are listed in Appendix [Cl To estimate 
the efficiency of these mechanisms, we compare the cool- 
ing time for each process, icooi = l/{d\nE/dt), to the 
Hubble time, tn = 1/H{z). Except for Compton scat- 
tering and photon-photon scattering, we approximate the 
cooling time by the mean free time as most of the energy 
is lost in the first interaction. If tn S> icoob the photons 
ionize the IGM, produce energetic electrons, or downscat- 
ter, very rapidly. Conversely, if tn ^ ^cooi, the universe 
is optically thin and most of the energy is lost through 
the redshifting of photons. The results of this comparison 
are shown in Fig. [T]and Fig. [21 At the relevant redshifts 
for hydrogen recombination, z ~ 700 — 1200, while the 
universe is not transparent at the relevant energies, it is 
also not sufficiently opaque that we can ignore redshift 
entirely. 

For photons with energies below ~ 10'^ eV and above 



FIG. 2: A comparison of the photon cooling time (from all 
processes) to the Hubble time over the entire redshift range of 
interest. The plot assumes a He mass fraction of 1/4, with a 
baryon density of 2.57 x 10~^ amu / cm'^ today, and the stan- 
dard ionization history and fiducial cosmology. The dashed 
line corresponds to fcooi = tn- There is a discrepancy between 
this figure and Fig. 2 in the originally published version of 
[37j : the authors of that paper have advised us that upon 
revising their calculation, their results now agree with ours. 



~ 10^^ eV, in the redshift range of interest, the domi- 
nant processes (photoionization and pair production on 
the CMB, respectively) take place on timescales much 
faster than the Hubble time. The lowest-energy photons 
deposit their energy into the IGM by photoionization, 
while the highest-energy photons rapidly pair produce 
or downscatter on the CMB. Photon-photon scattering 
is a "photon splitting" process that yields an approxi- 
mately fiat photon spectrum (up to the energy of the 
initial photon), whereas pair production produces an ap- 
proximately flat spectrum of high-energy electrons and 
positrons which rapidly inverse Compton scatter to pro- 
duce a softer photon spectrum. 

Photons lying in the broad 10'^ — 10^^ eV range scat- 
ter or pair-produce on timescales within a few orders of 
magnitude of the Hubble time, while slowly redshifting 
away their energy. With decreasing redshift, all the en- 
ergy loss processes decrease in efficiency relative to the 
Hubble time, as shown in Fig. ^ 

There is a "transparency window" at ^ 10^ — 10^° eV 
at z — 1000, where the cooling time of the dominant en- 
ergy loss processes is close to the Hubble time. The ratio 
tH/tcooi oc (1 -I- z)^/^ ((X (1 -I- z)^/^ for photon-photon 
scattering), while the photon energy redshifts as {1 + z): 
photons injected into a transparency window can there- 
fore remain in the optically thin regime, and contribute 
to the diffuse photon background today. Below this en- 
ergy range, Compton scattering rapidly depletes the pho- 
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ton spectrum, but becomes inefficient as an energy-loss 
process at lower energies where Compton scattering be- 
comes purely elastic. The result is a second, narrower 
"transparency window" below the Compton bump. 



C. Beyond the "on-the-spot" approximation 

Previous analyses of the effect of DM annihila- 
tion/decay on the ionization history of the universe 
d, [33, m, 113] have employed the "on-the-spot" approx- 
imation, where the energy from DM annihilation/decay 
is assumed to be instantaneously deposited in the IGM 
with some efficiency / ~ 0.01 — 1. The effects of Sommer- 
feld enhancement can be determined within these models 
simply by multiplying the DM annihilation cross section 
by the appropriate saturated enhancement factor. 

However, the on-the-spot approximation is not neces- 
sarily well justified at redshift ~ 1000. The energy loss 
processes for electrons, very low-energy photons and very 
high-energy photons are all extremely rapid compared to 
the Hubble time, but over much of the relevant energy 
range the timescale for photon energy loss is of the same 
order as the Hubble time (Fig. [2]). If a significant frac- 
tion of the energy from annihilation is injected at this 
energy scale or higher, assuming either that photons in 
the semi-transparent regime promptly deposit their en- 
ergy or that they never deposit their energy may give a 
poor approximation to the actual deposition history. 

Instead, we evolve the spectra of the photons from DM 
annihilation with the various energy loss processes de- 
scribed above, taking redshifting into account, and track 
the energy deposited as a function of redshift. However, 
at very high and low photon energies, the energy loss pro- 
cesses are very rapid compared to the Hubble time. Sta- 
ble numerical evolution of the photon spectrum requires 
that we integrate out these rapid processes (or take ex- 
tremely small timesteps, but this is not practical when 
the timescale for the fastest cooling mechanisms may be 
12 orders of magnitude shorter than the Hubble time). 

At very low energies we may simply assume prompt 
deposition with 100% efficiency, but resolving the high- 
energy processes is more complicated. Although a 100 
GeV photon will scatter or pair-produce very quickly rel- 
ative to the Hubble time, this does not imply that all 
its energy is deposited on that timescale: photon-photon 
scattering on the CMB will split the original photon's 
energy repeatedly until it enters the transparency win- 
dow, and inverse Compton scattering on the electrons 
and positrons resulting from pair production can also in- 
ject high-energy photons into the transparency window. 



D. Resolving rapid cooling processes for high 
energy photons 

To resolve the effect of high-energy photon cooling pro- 
cesses, we begin by ignoring redshifting entirely (since by 



construction, these processes occur on timescales several 
orders of magnitude shorter than a Hubble time). Our 
goal is to compute the total energy deposited to the IGM, 
and the resulting spectrum of lower-energy photons, once 
no photons remain at the highest energies, due to pair 
production and downscattering (which occurs quickly rel- 
ative to a Hubble time). 

We compare the rates of the various processes to de- 
termine the fraction of the initial photon's energy that 
is directly injected into the IGM, and the spectrum of 
lower-energy photons that is created, in the photon's first 
interaction. For high energy photons, the former occurs 
primarily by pair production followed by excitation, ion- 
ization or heating by the newly produced electrons or 
positrons. The latter occurs by (1) photon-photon scat- 
tering, and (2) pair production followed by inverse Comp- 
ton scattering and positronium annihilation of the e~^e~ 
pairs. As described previously, we integrate out the elec- 
tron cooling processes to determine the scattered photon 
spectrum and deposited energy resulting from pair pro- 
duction. 

Part of the resulting photon spectrum may still lie in 
the high energy region where the energy loss processes 
are too rapid (compared to tn) to be easily included 
when calculating the evolution of the spectrum over cos- 
mological timescales. Thus we need to apply the same 
procedure to the high energy bins in the new spectrum, 
determining the photon spectrum produced by downscat- 
tering and pair production from those bins. Iterating 
this process is equivalent to taking into account multiple 
interactions for the high-energy photons inside a single 
timestep, and the resulting infinite series can be written 
in terms of the sum of a geometric series of square ma- 
trices (of size equal to the number of high energy photon 
bins). The series can then be resummed analytically by 
a matrix inversion, which is computationally tractable so 
long as the number of high energy bins is not too large. 



E. Time evolution of the photon spectrum 

We employ a simple first-order integration scheme for 
the time evolution of the photon spectrum and energy 
deposited in the IGM. We initialize the photon spectrum 
as zero everywhere (we do not track the CMB photons), 
so there is an initial transient behavior due to the neglect 
of photons from earlier DM annihilations. However, at 
z > 2000 the deposition efficiency is excellent and the 
on-the-spot approximation is quite accurate, so the tran- 
sient behavior dies away rapidly provided the initial red- 
shift is sufficiently high. We choose an initial redshift of 
^^init ~ 4000. We employ the standard ionization his tory 
computed by the publicly available code RECFAST [50], 
and assume perturbations to the ionization history due 
to the extra injected energy to be small (since large per- 
turbations to the ionization history would violate existing 
constraints). 

At each timestep, we divide the photon spectrum into 
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three regions: low energy photons where the (inelastic) 
interaction rate is short compared to H{z) ("fast" defined 
as a total rate of more than one interaction per photon 
per timestep, for timcsteps typically of order ^ ijj/lOOO), 
high energy photons where the interaction rate is fast 
compared to H{z), and an intermediate region where the 
interaction rates are no more than a few orders of mag- 
nitude greater than H{z) and so the timestep is small 
enough to resolve them. We take a standard timestep of 
(iln(l + z) = 10"'^. This division is redshift-dependent 
and so must be performed at each timestep. The low 
energy photons are assumed to deposit all their energy 
within the timestep, and the high energy photon pro- 
cesses are integrated out as described above. The part 
of the photon spectrum in the intermediate energy range 
is evolved by the photon cooling mechanisms listed in 
Appendix [Cl and by redshifting. At the end of each 
timestep, the photon spectrum is updated and new pho- 
tons are injected from DM annihilation (both from direct 
production, and from ICS and annihilation of electrons 
and positrons). The energy deposited to the IGM at each 
step is recorded. 

The evolution of the photon spectrum is shown for 
two sample models in Fig. O The effect of the semi- 
transparent windows discussed previously is clear, with 
large peaks in the spectra at ~ 10^ and ~ 10^ — 10^" eV. 
The gap between the peaks is due to Compton downscat- 
tering rapidly depleting the photon spectrum in this en- 
ergy range, while above and below the semi-transparent 
windows, pair production on the CMB and photoioniza- 
tion, respectively, dominate. The edge from positron an- 
nihilation is visible at 511 keV. 
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III. ENERGY DEPOSITED TO THE IGM 



A. The efficiency factor f{z) 



FIG. 3: The photon spectrum as a function of energy at sev- 
eral redshifts for A/dm = 1000 GeV (top) and Mdm = 10 GeV 
{bottom), for xX ~* 4>4' followed by ^ e^e~ , with — 1 
GeV. 



The power per hydrogen atom injected by DM annihi- 
lation is frequently written in the form, 

CDM = 2/Mdm (1 + zf (1) 

where riDM.o and nH,o are the WIMP and H number den- 
sities at z = 0, respectively. In the on-the-spot approxi- 
mation, / is just the efficiency with which the WIMP rest 
mass energy liberated by annihilation is injected into the 
IGM ^22,]. Note that / is defined in terms of the WIMP 
mass, not the total energy of the electrons produced by 
annihilation: a large branching ratio to neutrinos, for 
example, results in a smaller value for /. 

In our previous paper on this topic [1], we made the 
simplifying assumption that / and (aAv) were indepen- 
dent of redshift, an a ppr oach that has also been adopted 
by other authors (sol. |43|. The Sommerfeld enhancement 
can cause (aAv) to vary with redshift, but as discussed 



previously, in the redshift range relevant to this prob- 
lem we can reasonably assume that the enhancement is 
saturated (although our numerical code could trivially 
accommodate a time-dependent (aAv)), 

(aAv) = S,na^{aAv)fo ■ (2) 

Here {aAv)io is the usual thermal relic freeze-out cross 
section of ~ 3 x 10^^^ cm'^/s, and Sma.x is the saturated 
Sommerfeld enhancement. 

Our detailed numerical computation of the energy de- 
position allows us to go beyond assuming a constant /, 
for the models of interest. Because of the changing trans- 
parency window (Fig. ^ the rate at which the photons' 
energy is absorbed by the IGM varies with z and with 
WIMP model, even in the on-the-spot approximation. 
Where the on-the-spot approximation breaks down, the 
delayed absorption of annihilation energy injected at ear- 
lier times can also alter the effective f{z) profile. To 
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some degree, these effects may cancel each other out: 
the universe becomes more transparent at low redshifts 
and this reduces deposition efficiency, but there are more 
photons present than in the on-the-spot approximation, 
due to photons injected at higher redshift that have not 
yet completely cooled. 

In the absence of the on-the-spot approximation, the 
physical meaning of / as an efficiency factor is not so 
clear, but a slight extension of Eq. [T]is still a very useful 
parametrization. We write. 



eDM{z) = 2f{z) 



S„ 



DM 



/lOO GeV\ , 

2/(^)^,„ax — TT (1 + zy 

x2 X 10"2''cV/s/H. 



(4) 



It is this energy injection that the CMB data constrain. 
Note that the injected energy is inversely proportional 
to the particle mass; more massive particles inject less 
energy into the IGM. 



(aAv) = 2 X 10~^^ cm'^/s has been proposed as an ex- 
planation for the positron excess observed by PAMELA 
[55|. In this model the neutralino annihilates predom- 
inantly to W bosons, and f{z) and the impact of DM 
annihilation on the CMB can therefore be immediately 
computed from our results. We calculate the spectra of 
e~^e~ and photons resulting from these annihilation chan- 
nels (including final state radiation) using PYTHIA. Fig. 
|4] displays f{z) as a function of redshift for the various 
WIMP masses and annihilation channels. 

In Appendix[A]we present accurate fitting functions for 
f{z) over the redshift range z — 300 — 1200, for all the 
annihilation channels under consideration. Our choice of 
redshift range is dictated by the impact on the visibil- 
ity function: above z ~ 1200 or below z ~ 500, changes 
to f(z) have minimal impact on the last scattering sur- 
face (see for example Fig. 3 in Less precise fits are 
given for higher and lower redshifts. Note that these fit- 
ting functions do not apply to models of decaying dark 
matter: while our numerical approach can be applied to 
decaying WIMP models with only trivial changes, we de- 
fer such an analysis to future work. 



B. DM annihilation channels 

As discussed previously, recent cosmic-ray anomalies 
have motivated models of WIMP annihilation to leptons 
with a large cross section. We compute f{z) for a WIMP 
annihilating to lepton pairs and charged pions, both di- 
rectly and via a new GeV-scale state (annihilation chan- 
nels of the latter type are denoted "XDM" ) . As a bench- 
mark, the mass of the new light state is taken to be 1 
GeV for electron, muon and pion final states, and 4 GeV 
for taus: however, because of the large mass hierarchy be- 
tween the WIMP and the GeV-scale state, the spectrum 
of the SM annihilation products is nearly independent 
of this choice of parameter. We investigate a range of 
WIMP masses for these annihilation channels, including 
the mass/channel combinations fitted to PAMELA and 
ATIC data in ;20"j. We also compute f{z) for WIMPs an- 
nihilating to pairs of SM particles, for benchmark WIMP 
masses of 200 GeV and 1000 GeV. 

We note that at sufficiently high Mdm, the photons 
from annihilation are injected at energies where pair pro- 
duction on the CMB is extremely rapid, and so the spec- 
trum of produced photons is entirely determined by the 
pair production -|- ICS cascade. Consequently, around 
the redshift of last scattering, /(z) becomes essentially 
independent of the WIMP mass at Mdm^I TeV. 

Taking linear combinations of these channels allows 
f{z) to be computed accurately for a wide range of DM 
models. In particular, models where the DM does not 
originate in thermal equilibrium can have large annihi- 
lation cross sections, and so can be strongly constrained 
by the CMB once f{z) is known precisely. The indi- 
rect detection prospects of such nonthermal neutralino 
models have been investigated in [5l|, [H, HI, H^l- Re- 
cently, a 200 GeV wino with annihilation cross section 



IV. CMB CONSTRAINTS 

CMB constraints on the energy injection from DM an- 
nihilation, assuming a constant /(z), have been previ- 
ously calculated in ^ and updated in Using the 
approximate values for /(z) given in TableU we can now 
relate these constraints to specific DM models. A full 
analysis of the CMB constraints would take into account 
the variation of /(z) with redshift, but since /(z) is slowly 
varying in the redshift range of interest, we employ the 
approximation of a constant / to estimate the degree to 
which the CMB can constrain models of interest, and 
defer a more complete analysis to future work. 

To a large extent, the effect of DM annihilation on the 
CMB depends only on the average of /(z) in the redshift 
range of interest z ~ 800 — 1000, not the details of /(z). 
We employed the public codes RECFAST and CAMB [Sg 
to derive the effect of DM annihilation on the TT, TE 
and EE angular power spectra, for several different /(z) 
profiles (with the widest possible differences in mass and 
annihilation channel) , normalizing the annihilation cross 
sections so that the deposited power e(z) would be the 
same in all cases if /(z) was replaced by its mean value 
between z = 800—1000. Fig. [5] demonstrates the results: 
even in an extreme case where the cross section is so large 
as to be strongly excluded by current constraints, the 
differences in the spectra due to differing /(z) profiles 
are small. 



A. Models fitting cosmic-ray excesses 

We focus here on models which fit the cosmic-ray ex- 
cesses measured by PAMELA, and in the case of higher- 
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FIG. 4: The "deposited power fraction" f(z) is the ratio of the power deposited in the gas (in the form of ionizations, excitations, 
and heating) to the mass energy hberated by WIMP annihilations. For electron channels, f(z) ~ 1 at high z, but other channels 
lose some fraction of their power to neutrinos, protons and neutrons. Upper left panel: direct annihilation to SM leptons. Upper 
right panel: direct annihilation to non-leptonic SM states ("light quarks" corresponds to 50 % annihilation to u quarks, 50 % 
to d quarks). Lower left panel: XDM-type models with annihilation through an intermediate 1 GeV state to electrons and 
muons. Lower right panel: XDM-type models with annihilation through an intermediate 1 GeV state to charged pions, and 
through an intermediate 4 GeV state to taus. The legend indicates the annihilation channel and the WIMP mass. The kink 
around z = 1700 is an artifact of an approximation made in RECFAST and has no impact on our results. 



mass WIMPs, also ATIC or Fermi. Boost factors and 
WIMP masses are taken from [g^l for the leptonic and 
XDM channels, and from for annihilation to W 
bosons. Fig. [6]displays the WMAP5 constraints on these 
models, and the region of parameter space that will be 
probed by Planck. 

In all cases, the models which fit the cosmic-ray ex- 



cesses are close to being ruled out by WMAP5, at 95 % 
confidence. The tension is greater for models which fit 
the ATIC excess, where the boost factors given in [2^ 
are already excluded. However, this result does not rule 
out these DM models as explanations for the ATIC ex- 
cess, due to astrophysical uncertainties in the required 
boost factor. For example, the local DM density is only 
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FIG. 5: CMB power spectra for three different DM annihilation models, with power injection normalized to that of a 1 GeV 
WIMP with thermal relic cross section and f — 1, compared to a baseline model with no DM annihilation. The models give 
similar results for the TT (left), TE {middle), and EE (right) power spectra. This suggests that the CMB is sensitive to only 
one parameter, the average power injected around recombination. All curves employ the WMAP5 fiducial cosmology: the 
effects of DM annihilation can be compensated to a large degree by adjusting Us and as 



known to within a factor of ^ 2 (which is then squared 
to determine the annihilation rate) , and density enhance- 
ments from local substructure could also contribute an 
0(1) boost to the cosmic-ray flux. The excess measured 
by Fermi requires generically smaller boost factors than 
ATIC, by a factor of ~ 2 — 3: such models are not ruled 
out by WMAP5 even without taking into account astro- 
physical uncertainties, but will be constrained by Planck. 

The degree of uniformity between the models should 
not be surprising, despite the wide range of masses and 
boost factors. The variations in f{z) between different 
channels arise in large part from the energy carried away 
by annihilation products other than photons and elec- 
trons - but these annihilation products also do not con- 
tribute to the cosmic-ray excesses measured at ATIC and 
PAMELA. The cosmic-ray excesses are more sensitive 
measures of the high-energy spectrum of the annihilation 
products than the CMB, whereas the CMB is sensitive to 
soft photons and electrons which may be absorbed into 
the background in cosmic-ray measurements, but to a 
first approximation both measurements are simply prob- 
ing the total power in electrons (at least when the power 
in photons produced by annihilation is small). 



B. Implications for Sommerfeld-enhanced DM 
annihilation 



and multiple light force carriers (e.g. [23|), but in this 
work we will consider only the simplest case. 

If the dark matter particle couples to a scalar media- 
tor (f) with coupling strength A, then the enhancement is 
solely determined by the dimensionless parameters. 



(v/c) 



a 



£0 



aMuu' 



(5) 



where a = A^/47r. In the limit where the (p mass goes to 
zero (e^ —^ 0), the enhancement to the annihilation cross 
section - denoted S ~ can be determined analytically, and 
S ^ Tr/sy at low velocities. For nonzero e^, there are two 
important qualitative differences. The first is that the 
Sommerfeld enhancement saturates at low velocity-the 
attractive force has a finite range, and this limits how 
large the enhancement can become. Once the deBroglie 
wavelength of the particle (Mdmw)""'^ exceeds the range 
of the interaction m^^, or equivalently once e„ drops 
beneath 6(4, the Sommerfeld enhancement saturates at 
~ ^ [2^ . The second effect is that for specific values 
of 60, resonances occur where the enhancement scales as 
~ instead of ^ l/^v, potentially increasing the en- 
hancement factor by several orders of magnitude. In the 
resonant case the velocity at which the enhancement sat- 
urates is also smaller than in the non-resonant case (for 
the same value of e^,). 



As described in the Introduction, the CMB has the po- 
tential to act as an especially sensitive probe of DM mod- 
els with Sommerfeld-enhanced annihilation. The sim- 
plest example of the Sommerfeld enhancement with a 
massive mediator is the case of WIMPs interacting via a 
Yukawa potential. More complicated models can contain 
small mass splittings among the dark sector particles. 



1. Saturation of the enhancement 

At first glance it might appear that our calculation 
would not apply to Sommerfeld-enhanced models, due 
to the variation of the enhancement with velocity, since 
we have assumed a constant (cr^f) with respect to z. 
However, for models which are not already ruled out 



10 



Ruled out by WMAP5 



1 1 



Planck 
forecast 



13 



CVL 



1 XDM 2500 GeV, BF = 2300 

2 n> 1 500 GeV. BF = 1 1 DO 

3 XDM 2500 GeV, BF = 1 000 

4 XDM e-e 1000 GeV. BF = 300 

5 XDM 4:4:1 1 000 GeV, BF = 420 
B e*e 700 GeV, BF = 220 

7 15D0GeV. BF = 56D 

8 XDM 1:1:2 1500 GeV, BF = 400 

9 XDM |i*ii 400 GeV, BF = 1 1 
1011*11 250 GeV. BF = 81 

1 1 W-W 200 GeV. BF = 66 
12XDMe-e 15DGeV.BF = 16 
13 e*e 100 GeV, BF= 10 



100 

DM Moss [GeV] 



FIG. 6: Constraints on the annihilation cross-section {aAv) 
the efficiency factor /. The dark blue area is excluded by 
WMAP5 data at 95% confidence, whereas the lighter blue 
area shows the region of parameter space that will be probed 
by Planck. The cyan area is the zone that can ultimately be 
explored by a cosmic variance limited experiment with angu- 
lar resolution comparable to Planck. Constraints are taken 
from [42] (Fig. 4). The data points indicate the positions of 
models which fit the observed cosmic-ray excesses, as fitted in 
[13, [H. Squares: PAMELA only. Diamonds: PAMELA and 
Fermi. Crosses: PAMELA and ATIC. Error bars indicate the 
factor-of-4 uncertainty in the required boost factor due to un- 
certainties in the local dark matter density (any substructure 
contributions are not taken into account). For models labeled 
by "XDM" followed by a ratio, the annihilation is through an 
XDM intermediate light state to electrons, muons and pious 
in the given ratio (e.g. "XDM 4:4:1" corresponds to 4:4:1 
annihilation to e^e~, fi^ fi~ and tt^tt"). 



by WMAP5 constraints, either the enhancement must 
be saturated over the redshift range in question (z ~ 
100 — 4000), or a or f{z) must be extremely small - in 
which case the model could not explain the cosmic-ray 
anomalies described in the Introduction. For the models 
of greatest interest, the enhancement 5* thus provides a 
constant boost factor to the annihilation cross section at 
z 1000, and our constraints apply directly. 

At redshift z, the CMB temperature is ^ 2.35 x 
10^'*(1 + z) eV. This places an upper bound on the tem- 
perature of the DM: however, after kinetic decoupling 
the DM temperature evolves adiabatically as T cx z^, 
and thus the WIMPs can be much colder than the pho- 
ton temperature. p3| suggests u/c ~ 10~^ at z ~ 1000 
for a 100 GeV WIMP. 

If the enhancement is still unsaturated at such low ve- 
locities, then the force carrier must be extremely light 
compared to the WIMP mass. For the models recently 
proposed in the literature [111, [H, [H, , the enhance- 
ment has always saturated by this point as the force carri- 
ers are much heavier than 10~^Mdm- Other constraints 
on models with very low-mass mediators also exist: as 



one example, a 1/w enhancement which saturates at too 
low a velocity can also cause runaway annihilations in 
the first DM halos at the onset of structure formation 
(58| . Furthermore, as shown in Fig. [51 models which fit 
the recently observed cosmic-ray anomalies are already 
close to being ruled out by WMAP5. If the Sommer- 
feld enhancement in such models has not saturated by 
(v/c) ~ 10~^, this implies an effective cross section at re- 
combination ~ 4 — 5 orders of magnitude higher than in 
the present-day Galactic halo. Such models are therefore 
strongly excluded by WMAP5. Similarly, if the WIMP 
annihilates to the same particle which mediates the Som- 
mcrfeld enhancement, then in order for the enhancement 
to evade the constraints in Fig. [SI the coupling a between 
the WIMP and the force carrier must be extremely small 
- reducing the annihilation cross section at freeze-out to 
unacceptable levels for a thermal relic. Thus for a broad 
range of well motivated models, it is self-consistent to as- 
sume that the Sommerfeld enhancement is saturated for 
the redshift range of interest (z ~ 100 — 4000). 

We can write the 95 % confidence limits from WMAP5 
in terms of constraints on the total cross section, 



(cAw) saturated < 



3.6 X 10-24cni3/s fMuMC- 



f 



ITeV 



(6) 



or as constraints on the maximum saturated enhance- 
ment, relative to the thermal relic cross section (cr^u) — 
3 X 10-26 cm^/s. 



5„ 



< 



120 fMBMC^ 
f 



ITeV 



(7) 



In both cases values of / for the different channels are 
given in Table [11 

These results directly limit the maximum boost fac- 
tor possible from substructure, in Sommerfeld-enhanced 
models. There has recently been considerable interest 
in possible annihilation signals from dark matter sub- 
halos, where the DM velocity dispersion is reduced and 
the Sommerfeld-enhanced cross section is boosted (e.g. 
[59I . l60l . HH, 1121). However, the saturated cross section 
cannot be much larger than that required to fit the cos- 
mic ray anomalies, so for models which fit the cosmic ray 
anomalies, the lower velocity dispersion in subhalos will 
not result in a higher annihilation cross section. 



2. Sommerfeld-enhanced models fitting cosmic ray excesses 

In Sommerfeld-enhanced models which produce the ob- 
served excesses in e'^e^ cosmic rays, the saturation of 
the enhancement is even more constrained than in the 
general case. Since the cross sections required to fit 
the cosmic ray anomalies are already nearly excluded by 
WMAP5, as shown in Fig. [6l the enhancement must al- 
ready be close to saturation at ~ 150 km/s (5 x IQ-^'c), 
the estimated local WIMP velocity dispersion. Astro- 
physical uncertainties - in the propagation of cosmic rays. 
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the local dark matter density and the local velocity dis- 
tribution - can weaken this constraint by a factor of a 
few. 

In the case of the simplest Sommerfeld enhancement 
scenario, where the light force carrier generates a Yukawa 
potential between the WIMPs, this saturation require- 
ment disfavors models with <C 10~'*Mdm- Stronger 
constraints on the force carrier mass may be possible in 
the case of resonant enhancement, where e„ < at sat- 
uration. 



3. Annihilation through the force carrier 

A related but independent constraint on the force car- 
rier mass occurs for models where the WIMP annihilates 
to the same particle that mediates the Sommerfeld en- 
hancement, provided the WIMP is a thermal relic and 
constitutes 100 % of the dark matter [t^. In this case 
the WIMP mass and jts coupling to the light force carrier 
satisfy the relation |24| , 



a : 



10" 



A/dmc^ 



(8) 



270 GeV 

since the freeze-out cross section for WIMP annihilation 



to the light force carrier is (aAv) ~ Tra^/Mpj^. But in 



the case of non-resonant Sommerfeld enhancement with 
behavior similar to the simple Yukawa case, the saturated 
enhancement is given by, 

(Mdmc2)2 



S„ 



10- 



DM 



(270 GeV)m0c2 



(9) 



Applying our constraints on S'max in the form given in 
Eq. [3 it follows that. 



Mdm' 



>io- 



0.3 



(10) 



In this case the constraint does not depend on the 
cosmic-ray data and thus there are no uncertainties from 
present-day Galactic astrophysics. However, 0(1) model- 
dependent factors in the saturated enhancement and the 
annihilation cross section at freeze-out may modify Eq. 
1101 and in the case of resonant enhancement, the con- 
straint can become significantly stronger (by a factor 
Smax^<p)' these factors can be computed given a specific 
DM model, and such a calculation is necessary to make 
the limit precise. As previously, though, force carriers 
with masses ^ 10^*7\fDM are disfavored. For WIMPs 
with masses of order a few TeV, as suggested by fits to 
the Fermi and H.E.S.S. data, this constraint may be suf- 
ficient to rule out force carriers with masses below the 
muon threshold, in a wide range of models. 



V. CONCLUSIONS 

The CMB provides a robust constraint on models of 
dark matter annihilation, which is particularly relevant 



for models that annihilate rapidly enough today to pro- 
duce the PAMELA, Fermi and/or ATIC signals. These 
signals all require an annihilation cross section larger 
than the thermal relic cross section, implying either non- 
thermal production or some other mechanism, such as 
Sommerfeld enhancement. In either case, the implied 
cross section at z « 1000 is as large as today, or possibly 
larger. 

Extending previous work d, [s^, IH, El] , we have com- 
puted the efficiency of annihilation power deposition in 
the IGM for a number of models, some tuned to fit the 
recently measured cosmic-ray anomalies, and some with 
more general annihilation channels. We have improved 
previous calculations by finding and correcting some er- 
rors in the photon cooling processes, and by treating the 
redshifting and downscattering of photons in the semi- 
transparent windows in detail. Numerical convergence 
on this calculation indicates it is correct at the ^ 1% 
level, an accuracy far exceeding that of any possible con- 
straints on local DM annihilation, as they depend on e.g. 
the local DM density squared, which is uncertain. 

Given the state of the current data, it is adequate 
to approximate the efficiency of energy deposition, f[z), 
with an average over the redshift range of interest. Fu- 
ture improvements in the data, however, may demand a 
more detailed calculation. Therefore we provide fitting 
functions for 40 models of interest in Appendix \K\ 

WMAP5 already constrains Sommcrfcld-enhanced 
models that have been proposed to explain the recently 
observed cosmic-ray excesses, requiring that the satu- 
rated annihilation cross section does not greatly exceed 
the cross section in the neighborhood of the Earth. This 
result argues against large DM annihilation signals from 
substructure, in the present-day Galaxy. In the sim- 
plest case where the force carrier is a scalar generating 
a Yukawa potential, and either the model produces the 
observed cosmic-ray signals or the DM is a thermal relic 
annihilating to the same particle that mediates the Som- 
merfeld enhancement, masses for the force carrier much 
less than 10~*Mdm are disfavored. 

For a local DM density of 0.3 GeV/cm'^, standard as- 
sumptions for cosmic-ray propagation [20| , and no contri- 
bution to the annihilation flux from substructure, there 
is mild tension (at the factor of 2 level) between the 
WMAP5 95 % confidence limits and models fitting the 
excess observed by ATIC. However, due to 0(1) uncer- 
tainties in the local density and substructure contribu- 
tion, these explanations for the ATIC excess cannot be 
excluded by WMAP5. Models fitting the excesses mea- 
sured by PAMELA and Fermi are consistent with the 
WMAP5 constraints, for standard astrophysical assump- 
tions. However, all of these models will be rigorously 
tested by the recently launched Planck experiment. 
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APPENDIX A: FITTING FUNCTIONS FOR f{z) 



Around the redshift of last scattering, f(z) is slowly varying with respect to redshift; the approximation of a constant 
f{z) is not too far from the truth. As described in Section ITVl we average /(z) over the redshift range z = 800 — 1000 
to obtain an effective constant efficiency /, for a range of masses and annihilation channels. For applications where 
a more accurate approximation to f{z) is required, we fit f{z) to a seven-parameter function over the redshift range 
z = 300 - 1200, 

/,.) = f (1 + =)• f ri±iV + f V' V =p ( ' . ) (Ai) 



l + ((l + z)/1100)'' 



These fits are accurate to within 1% between z = 300 — 1200 for all channels. These fits remain accurate to < 5% 
between z — 170 and z — 1470, but outside this range they may perform very poorly. 
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27 





31 





0809 


2 


5357 


4 


7587 





0457 





3035 


-0 


3322 





5392 





0179 


-13 


3422 


321.8945 


followed by 


400 





25 





31 





0741 


11 


3064 


3 


3949 





0402 





2937 


-0 


2579 





5965 





0505 


-10 


3800 


774.7615 




1000 





25 





31 





0617 


12 


5195 


3 


1133 





0418 





2925 


-0 


3294 





7487 





0541 


-10 


6936 


939.3080 




2500 





24 





31 





0556 


13 


0389 


2 


9343 





0522 





2926 


-0 


6537 





5413 





0566 


-10 


5987 


952.4342 


XDM taus 


200 





22 





27 





0604 


6 


6206 


3 


6373 





0333 





2861 


-0 


0610 


1 


0364 





0548 


-8 


7336 


638.6944 


XX —' 4>4-, 4> -r + T~ 


1000 





22 





27 





0534 


11 


2208 


3 


1869 





0424 





2841 


-0 


4351 





6734 





0542 


-10 


5137 


911.3169 


XDM pions 


100 





22 





25 





0607 


1 


4685 


5 


0403 





0394 





2881 


-0 


2700 





5445 





0137 


-12 


6965 


304.5202 


XX 0<A 


200 





21 





25 





0674 


6 


0060 


4 


1253 





0353 





2825 


-0 


1722 





7910 





0323 


-13 


6145 


477.7644 


followed by 


1000 





20 





25 





0515 


12 


3319 


3 


1745 





0382 





2762 


-0 


3601 





6781 





0517 


-10 


8809 


1030.3075 




1500 





20 





25 





0481 


12 


6927 


3 


0715 





0428 





2760 


-0 


5297 





5865 





0547 


-10 


7564 


1026.1082 




2500 





20 





25 





0453 


12 


9871 


2 


9688 





0480 





2762 


-0 


6968 





5217 





0566 


-10 


6509 


1025.4334 


W bosons 


200 





29 





35 





1013 


19 


1565 


2 


9322 





0395 





3076 


-0 


0895 


1 


1093 





0377 


-13 


2287 


446.3091 


XX VU+ W 


300 





29 





35 





0906 


15 


7615 


3 


0067 





0388 





3053 


-0 


0855 


1 


0554 





0389 


-13 


1812 


528.0655 




1000 





28 





35 





0711 


10 


6406 


3 


1935 





0415 





3025 


-0 


2181 





8366 





0516 


-10 


0585 


782.1619 


Z bosons 


200 





28 





34 





0998 


20 


7336 


2 


8932 





0392 





3043 


-0 


1088 


1 


0375 





0359 


-13 


3227 


447.9354 


XX zz 


1000 





27 





33 





0689 


10 


6396 


3 


2027 





0407 





2988 


-0 


2263 





7934 





0514 


-9 


9893 


773.0394 


Higgs bosons 


200 





34 





40 





1313 


24 


2160 


2 


8491 





0479 





3205 


-0 


2349 





7599 





0297 


-13 


5576 


388.8721 


XX — * /i/j- 


1000 





32 





40 





0877 


10 


9585 


3 


1982 





0430 





3133 


-0 


1570 





8487 





0490 


-9 


8120 


616.1287 




200 





35 





41 





1244 


20 


6286 


2 


8789 





0467 





3217 


-0 


1873 





8494 





0345 


-13 


3583 


383.5586 


XX — * fafa 


1000 





33 





41 





0917 


11 


6611 


3 


1846 





0425 





3149 


-0 


1246 





9724 





0467 


-9 


8366 


635.3690 


Light quarks 


200 





34 





40 





1129 


18 


5995 


2 


9221 





0432 





3174 


-0 


1218 





9244 





0361 


-13 


1747 


430.2257 


XX ^ "O. (50 % each) 


1000 





32 





40 





0882 


12 


3648 


3 


1280 





0434 





3135 


-0 


1700 





9101 





0490 


-9 


8913 


674.5797 



TABLE I: Column 3: Mean values of f{z) averaged from z — 800 — 1000. Column 4: f{z = 2500), where / approaches the 
asymptotic high-z limit. Columns 5-7: Fit parameters for Eq. IA2I the resulting fit is accurate to within 5% for z = 10 — 170, 
up to the limits of this calculation. Columns 8-14: Fit parameters for Eq. lAll the resulting fit is accurate to within 1% for 
z = 170 - 1470. 



At very high redshifts, z 3> 1000, f{z) converges to asymptotic values determined solely by the annihilation channel, 
since in this limit all the energy injected in electrons and photons (as opposed to neutrinos, protons etc) is efficiently 
deposited. In the (relatively) low-redshift region, z = 10 — 170, f{z) can be fitted to within 5% for all channels 
considered by a simple three-parameter fit. 



/(z) =aexp((ln((l + z)/6)/c)' 



(A2) 



However, at z ^ 100, the onset of structure formation may invalidate this calculation. 

Table H] lists the averaged / employed in Section HVl the value of f{z) at z = 2500, and the fit parameters for Eqs 
IAmA21 for the masses and annihilation channels discussed in Subsection IIII Bl 
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APPENDIX B: ELECTRON COOLING MECHANISMS 
1. Inverse Compton scattering 



The cooling time for ICS is [461, 



where Tcmb — 2.725(1 + z)K is the mean CMB temperature at the relevant redshift, aji is the radiation constant, 
and (TT = 6.65246 x 10^^''' cm^ is the Thomson cross section. Comparing this to the Hubble time, one finds 

i / 1 I \ 5/2 1 



10-^ -TKKK -1^1 ' (B2) 



icooi V 1000 J \/rw? 

Thus the timescale for ICS is much shorter than the Hubble time. 

At the highest energies 10 GeV and higher, depending on the redshift) the Thomson cross section formula is 
no longer valid. Consequently we must use the general form for the photon spectrum resulting from ICS [46| . For 
convenience, all energies are written in units of the electron mass, WgC^ ~ 511 keV. The doubly differential spectrum 
is given by, 

'^^ ^(JTC^{2q\ogq+{l + 2q){\-q) + 0.h{\^q)(Tqf/{l + Tq))n{t), (B3) 



dE^dedt 4 eE'i 



q ■ 



1 



E, r(l - E.JE,) 



Here e is the energy (before upscattering) of the soft photon that the electron scatters on: in this case, a CMB photon. 
Ef, is the electron energy, i?^ is the energy of the upscattered gamma ray, and rt(e) describes the energy distribution 
of the soft photons per unit volume (i.e. n(e) = dNcMB / dedV) . This spectrum must then be integrated over the 
energy distribution of CMB photons to determine the total photon spectrum from ICS of a high-energy electron on 
the CMB. 

2. Ionization, Excitation and Collisional Heating 

We use the fits presented by Arnaud and Rothenflug [6^ for ionization of neutral hydrogen, neutral helium and 
singly ionized helium, 

aiE.) = lO-W^;^^ |a (l - i) + (l - i)% Clnu + Dlnu/u^ (B4) 

where / is the ionization potential, E^ is the kinetic energy of the electron and u = Ee/I- The coefficients A, B, C, 
and D are given by, 

• Hydrogen: A = 22.8, B = -12.0, C = 1.9, D = -22.6, with / = 13.6 eV. 

• Neutral helium: A ^ 17.8, B = -11.0, C = 7.0, D = -23.2, with / 24.6 eV. 

• Singly ionized helium: A = 14.4, B ^ -5.6, C = 1.9, D = -13.3, with / = 54.4 eV. 

For excitation of hydrogen and neutral helium we use the fitting functions given by Stone, Kim and Desclaux [63 |. 
For kinetic energies well above threshold, the cross sections have the form, 

^cH.cHe = ^ , y , ^ (A HE^iJR) + B + CR/E^,^) . (B5) 

J^kin + -C/bin + ^cxc 

Here R w 13.6 eV is the Rydberg energy, oq w 0.529 x 10~^ cm is the Bohr radius, E'bin is the binding energy of 
the electron to be excited, £^exc is the excitation energy, and i?kin is the kinetic energy of the incident electron. A, B 
and C are dimensionless parameters with values given by. 
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• Hydrogen: A = 0.5555, B = 0.2718, C = 0.0001, 

• Neutral helium: A = 0.1771, B = -0.0822, C = 0.0356. 

At high energies the result for hydrogen agrees well with that of ShuU and van Steenberg [l^l, 

2.75 X 10-15 ln(£;/13.6eV) , ,^r.^ 
aeH ^ ^cm^ . (B6) 

For excitation of singly ionized helium we follow Fisher et al [65| . 

9 

o'oHo+ — TTT- {oi ln.x + pinxlx + 7 + 5/x + i]/x^) , (B7) 
lox 

a = 3.22, /? = 0.357, 7 = 0.00157, (5 = 1.59, r] = 0.764. 

Collisional losses become important at lower energies, 100 eV < E <1 keV. The cross section for collisional heating 
is given by 

(Tee = (7.82 X 10-")(0.05//)(lnA)(£;e/eV)-2cm2, (B8) 

where / = l^E/E ~ 0.05 is chosen to simulate the discrete nature of Coulomb collisions. We take the Coulomb 
logarithm In A ^ 10. 

APPENDIX C: PHOTON COOLING MECHANISMS 

1. Pair production on the CMB 

The doubly differential electro n sp ectrum for pair production by a gamma ray (energy E^) encountering a soft 
photon (energy e) is given by [66l l67j | , 

dN 3 1 ( AAHt,{UEM-E.)/A) ^ ^ , 2{2eA-l)A^ ( l\ A^ \ 

-^TC—^-^\ -— — 8eA+—— —- [I - — ]-— — ^ n(e) (CI) 



dEededt 64e2£;M Ee{A - E^) Ee{A~ E^) \ eAjEl{A~EeY 




Here A = E^ + e. As in Eq. IB3[ all energies are in units of the electron mass, and n(e) describes the energy distribution 
of the soft photons per unit volume. This expression must be integrated over the CMB energy distribution to obtain 
the spectrum of electrons produced by pair production from a high energy gamma ray on the CMB. The positron 
spectrum is identical. 

2. Pair production on H and He 

We employ the high energy (Born approximation) cross sections for pair production on ionized H, free electrons, 
and singly ionized He, as described by Motz, Olsen and Koch [6^. As previously, we write all energies in units of the 
electron mass. In the high energy limit (complete screening). He"*" and H"'' can both be regarded as singly charged 
point charges, so share a single cross section for pair production, 

a^arl{^H2E,)-'^Y (C2) 



where a is the fine structure constant and tq = 2.8179 x 10 cm is the classical electron radius. If this equation is 
also used to describe pair production in the field of free electrons, then in the limit of full ionization (including double 
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ionization of helium), this expression agrees with that given in [4^ for a fully ionized medium. However, we slightly 
modify this formula as suggested by Joseph and Rohrlich for pair production off electrons, 

a^arl(^H2E,)~'-^]. (C3) 
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For pair production on neutral hydrogen and helium, we employ the cross sections given in [48| . 

In all cases we use the high energy Born approximation form for the spectrum of produced pairs (ssj . 

3. Photon-photon scattering 

Photon-photon scattering occurs when a gamma ray upscatters a CMB photon. It can be regarded as a photon 
"splitting" process, with each photon in the final state carrying away 1/2 the energy of the original gamma ray. 
We use the rate and spectrum given in [i^. 

The total rate for scattering of a gamma ray of energy on the CMB (already integrated over the CMB spectrum) 
is given by, 

R = 1.83 X IQ-^^ h^o^Tl^H„{^ + zfE^, (C6) 

where T2.7 is the present CMB temperature divided by 2.7 K, is measured in units of the electron mass, and 
/150 = 7?o/(50km/s/Mpc) ~ 1.4. The spectrum of resulting gamma rays is given by, 

«(.0/7)±fl-#+f#)T. ,C7) 



dEdt E^ \ E.y y E-^ 



4. Compton scattering 



Following [37|, we treat all electrons as free for Compton processes. The differential cross section for Compton 
scattering is just the usual Klein-Nishina cross section [70], given by, 

where Tc — 3.862 x 10^^^ cm is the Compton radius of the electron, E^ is the initial energy of the photon, E'^ is its 
final energy, and all energies are in units of the electron mass. 

5. Photoionization 



We include photoionization on neutral hydrogen, neutral helium and singly ionized helium. It is crucial to include 
photoionization on helium, as otherwise photoionization appears to turn off at rcdshifts greater than z ^ 1000 (i.e. 
hydrogen recombination), which does not reflect the true situation. 
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The photoionization cross sections for hydrogen and singly ionized heUum are known analytically, and are given by 

El, 

1 



thrcs J- 

297r2r| ^^^4 exp(-4??arctan(l/??)) 
" = ^^(^thrcs/i?) i_exp(-27rr;) ' ^^^"^ 

where tq — 2.818 x 10""'^^ cm is the classical electron radius, and -Ethrcs — 13.6 eV for H, 54.4 eV for Hc+. 

Svensson and Zdziarski [isj fitted the photoionization cross section for a mixture of 75% neutral hydrogen and 25% 
neutral helium (by mass) with a broken power law. We subtract the hydrogen contribution and use the result as the 
photoionization cross section for neutral helium, 

.(He) = -12.(H) + 5.1 X lO^W ( ^^^Z ' f <f '^"^ (CU) 

I (ssiv) ' ^ > 250eV. 
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